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Abstract. We extend earlier work on traveling waves in premixed flames in a grav- 
itationally stratified medium, subject to the Boussinesq approximation. For three- 
dimensional channels not aligned with the gravity direction and under the Dirichlet 
boundary conditions in the fluid velocity, it is shown that a non-planar traveling 
wave, corresponding to a non-zero reaction, exists, under an explicit condition relat- 
ing the geometry of the crossection of the channel to the magnitude of the Prandtl 
and Rayleigh numbers, or when the advection term in the flow equations is neglected. 



The Boussinesq-type system of reactive flows is a physical model in the descrip- 
tion of flame propagation in a gravitationally stratified medium [23]. It is given as 
the reaction-advection-diffusion equation for the reaction progress T (which can be 
interpreted as temperature), coupled to the fluid motion through the advection ve- 
locity, and the Navier-Stokes equations for the incompressible flow u driven by the 
temperature-dependent force term. After passing to non-dimensional variables [31 [19] , 
the Boussinesq system for flames takes the form: 



Here, z/ > is the Prandtl number, that is the ratio of the kinematic and thermal diffu- 
sivities (inverse proportional to the Reynolds number). The vector p = pg corresponds 
to the non-dimensional gravity g scaled by the Rayleigh number p > 0. The reaction 
rate is given by a nonnegative 'ignition type' Lipschitz function / of the temperature, 
this last one normalized to satisfy: < T < 1. The above model can be derived from 
a more complete system under the assumption that the Lewis number equals 1. 

We study the system (11.11) in an infinite cylinder D gM.^ with a smooth, connected 
crossection C M^. Recent numerical results, motivated by the astrophysical context 
[m [20] , suggest that the initial perturbation in T either quenches or develops a curved 
front, which eventually stabilizes and propagates as a traveling wave. On the other 
hand, existence of non-planar traveling waves for the single reaction-advection-diffusion 
equation in a prescribed flow has been a subject of active study in the last decade 
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p3t El [11]. For system (11. ip . existence of traveling waves has been considered under 
the no-stress or Dirichlet boundary conditions in u, in channels of various inclinations 
and dimensions H [3l [13 El [lOj . 

The main difference presents itself at the orientation of D with respect to g; when 
they are aligned there are no non-planar fronts at small Rayleigh numbers while 
in the other case a traveling front, necessarily non-planar, is expected to exist at any 
range of parameters. This has been rigorously proven: in [3] for n = 2 dimensional 
channels D and under no-stress boundary conditions, in [S] for n = 2 and the more 
physical no-slip conditions and in [TU] for the same boundary conditions and arbitrary 
dimension n, but for a simplified system (corresponding to the infinite Prandtl number 
u = oo) when the Navier-Stokes part of (II. ip is replaced by the Stokes system. 

The purpose of this paper is to remove this last assumption, for three dimensional 
channels. Namely, we will investigate the model supplied by the Navier-Stokes system. 
We assume that p is not parallel to the unbounded direction of D, which after an an 
elementary change of variables [3] amounts to studying: 

D = (— oo, oo) X Q = {{x,x); x G M, x e Q} 

and 

We will prove the existence of a traveling wave solution to (II. ip : T{x—ct, x), u{x—ct, x), 
with the speed c to be determined and under the boundary conditions: 

dT 

(1.2) — = and M = on dD, 

on 

where n is the unit normal to dD. Such a front satisfies: 

-cT^-AT + u-VT = f{T) 

(1.3) —cUx + du ■ Vm — uAu + Vp = Tp 

div M = 0. 

We set constant d to be or 1. For the simplified system, when the advection in u has 
been neglected and d = 0, the theorem below states existence of a non-planar traveling 
wave, for any crossection f2, Prandtl number u and Rayleigh number p. 

For the full system when d = 1, we need to assume the following relative thinness 
condition, involving z/, p, the area and the Poincare and the Poincare-Wirtinger 
constants Cp, Cpw of VL\ 

(1.4) v^-^|^]|l/2 |^|p1Cp^+ (^£|p-(0,x)p^ ^ j <1. 

This condition is essential in our analysis and it is not clear if the below existence result 
holds without it. Recall that Cp is determined by the thinness of f2, and hence (11.40 
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admits domains with large area which are sufficiently thin. Respectively, Cpw depends 
on the maximum of (inner) distances between points in Vt. On the other hand, the 
quantities relating to smoothness of dVL have no direct influence on fll.4p . 

The nonlinear Lipschitz continuous function / is assumed to be of ignition type: 

f{T) = on (-00, ^0] U [1, 00), f{T) > on (^o, 1) 

for some ignition temperature Oq G (0, 1). 
The following is our main result: 

Theorem 1.1. Assume that either d = or d = 1 and holds. Then there exist 

c>0,T e C'^^^{D) with VT G L'^{D), u e f] C'^^^{D), p G C^^iD) satisfying / fO) 
and together with: 

(1.5) lim ||m(x, Ollc^cn) = 1™ ||VT(x, ^lU-cn) = 0. 

X— >±oo x— >±oo 

Moreover T{D) C [0,1], maXa.>o,ygr2 T(x, y) = 6q, and there is a nonzero reaction: 

[ /(T)G (0,00). 
Jd 

The limits ofT satisfy: 

lim I |T(x, ■) I |L°o(m = 0, lim | \T{x, ■) - ^- 1 |L=°(n) = 

2;— >+oo x— >— 00 

for some: 6'_ G (0, 6*0] U {1}. 

The following sections are devoted to the proof of Theorem 11.11 In section 2 we 
formulate some auxiliary results, of an independent interest. In particular, we prove a 
weak version of Xie's conjecture [22] for the Stokes operator (established in [21] for the 
Laplacian). Based on results in [H], we then derive an a priori estimate valid in any 
channel D, whose cross-section Q fulffils the geometrical constraint (11. 4p . This allows 
us to obtain uniform bounds on the quantities involved in the fixed point argument 
(Theorem 14.21) in sections 3 and 4; in particular the bounds are independent of length 
of the compactified domains Ra = [—a, a] x f2. The set-up for the Leray-Schauder 
degree is different than in [5l |10]: we solve the flow equations in the full unbounded 
channel D, while the reaction equation is solved in Ra- Once the uniform bounds 
are established, we refer to [5], (TU] for further details of the proofs. In section 5 we 
improve a sufficient condition from [TO] for the left limit 6^ of the temperature proffie 
T obtained in Theorem II. II to be equal to 1. 

We remark that the a priori estimates we derive do not preclude the solutions (T, u) 
to have arbitrary large norms. Indeed, the main chain of estimates eventually leads 
to inequality (14.41) . whose right hand side has a linear growth in terms of the left 
hand side, and thanks to condition (II. 4p . the main bound Oil 1 1 W 1 1 7^00 does not restrict 
the magnitude of this quantity. Similar estimates are known also for solutions to the 
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Navier-Stokes equations for the 2d and cylindrical symmetric systems [HI UHl [12], [13] , 
and in presence of a special geometrical constraint on the domain fTIj [T3] . 

We will always calculate all numerical constants at the leading order terms explic- 
itly. By convention, the norms of a vector field u on D are given as: ||'u||l°°(_d) = 

(ELi ^ and \\u\\l2^d) = (EL h'Wl^iD)) ' ■ 

Acknowledgments. M.L. was partially supported by the NSF grant DMS-0707275 
and by the Center for Nonlinear Analysis (CNA) under the NSF grants 0405343 and 
0635983. P.B.M. has been supported by MNiSW grant No. N N201 268935 and by 
ECFP6 M. Curie ToK program SPADE2, MTKD-CT-2004-014508 and SPB-M. 



2. Auxiliary results 
In the sequel, we will need a uniform estimate for the supremum of the solution to 

1 /2 1 /2 

Stokes system in Z). The known proofs of the inequality IImUlo" < C'q||Vm||/2 ||^^'w||^2 , 
V being the Helmholtz projection, are based on the a-priori estimates in [1], which hold 
for smooth domains. Therefore the constant depends strongly on the boundary 
curvature, and becomes unbounded as tends to any domain with a reentrant corner. 
It has been conjectured by Xie [22] that actually = 1/ v^Svr. To our knowledge, this 
is still an open question. Below we prove its weaker version, sufficient to our purpose 
and involving lower order terms. 

Theorem 2.1. Let g G L'^{D). Then the solution u e H'^ H H^{D) to the Stokes 
system: 

(2.1) —uAu + 'Vp = g, div u = in D 

satisfies the bound: 

where constant Cq depends only on the crossection Q. 

Proof. 1. We first quote two results, whose combination will yield the proof. The first 
one is Xie's inequality [21j for the Laplace operator in a 3d domain. Namely, for any 
ueH'^r\ Hl{D) there holds: 

(2-2) \\u\\l^{d) < -^||An||J^(^(^JVu||J^(^.^.. 

The crucial information in the above estimate is that the constant 1 / \/2tx is good for 
all open subsets of M^. 
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The next result is a recent commutator estimate by Liu, Liu and Pego Recall 
first [12] that for any vector field u G L'^{D) there exists the unique decomposition 
u = Vu + Vg with dwiVu) = 0, and q solving in the sense of distributions: 

dq 

Aq = div u in D, — = on dD. 

on 

This Helmholtz projection satisfies: < ||m||l2(d)- In this setting, it has been 

proved in pjj that for every e > there exists C^^n > such that: 

(2.3) yu e n H^iD) J |(AP-PA)m|2 < Q + e^ J \Au\^ + C,,n J |V 

The proof in [JJJ, written for bounded domains, can be directly used also for the case 
of cylindrical domains D with smooth boundary (since the covering number for the 
partition of unity on dD is finite). 

2. Applying the Helmholtz decomposition to (12.11) we arrive at: —uVAu = Vg, 
which can be restated as: 

-Am = (VA - AV)u + ^Vg, 

since Vu = u. Using (12.31) we obtain: 

3 1 
||Am||l2(£,) < -||AM||i2(£,) + Cf^ll VM||i2(£,) + -||P5'||l2(d), 

which yields: 

4 

(2.4) \\Au\\l2(j^) < -\\g\\L^(D) + e'en Vn||i2(£,). 

Now combining (12.41) and (12.21) proves the result. ■ 

We will also need an extension result for divergence free vector fields. Define a 
compactified domain Ra = [—a, a] x Q. 

Theorem 2.2. For any a > and any e > there exists a linear continuous extension 
operator E : C^'°'{Ra) — > C^'°'{D), such that for every u G C^'°'{Ra) there holds: 

(i) {Eu)\R^ = u, 

(ii) if div u = in Ra, then div (Eu) =0 in D, 

(iii) \\Eu\\l--^d) < (1 + e)||M||Loo(R^). 

Proof. Given a vector field u G C^'"([— a, 0] x Q) we shall construct its extension u G 
C"^'°([— a, oo) X Q) such that (ii) holds together with: 

(2.5) II u\\L°°{[-a,oo)xn) < (1 + e) ||'if||L°°([-a,0)xn). 

This construction, being linear and continuous with respect to the C^'" norm, will be 
enough to establish the lemma. 
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Fix a large n > 0. For x G [— a,a/2n^] and x E fl, define the vector v{x,x) with 
components: 

u^{x,x) for X G [— a, 0] 

\iu^{0, x) + X2U^{—nx, x) + X^u^^—n^x, x) for x G [0, a/2n^], 



.1/ 



for i = 2, 3 : 



where: 



_ u\x,x) for X G [— a, 0] 

^ ' ' —n\2u'^{—nx,x)—n'^X-iU^{—n'^x,x) for x G [0, a/2n^], 

_{l + n){l+n^) 1 + ^2 1 + n 

Xl — 7, , A2 — ^7 -T-, A3 — — — -. 



Since we have: Yl'i=i Xi = 1, —nX2 — n'^X-s = 1 and 11^X2 + n^X^ = 1, it follows that 
V G C^'"([— a, a/2n2] x Q). Also, by an explicit calculation, we see that div -u = 
implies div v = 0. 

Let now G C°°(M, [0, 1]) be a non-increasing cut-off function such that (f){x) = 1 
for X < and (f){x) = for x > a/Sn^. Define: 

u{x, x) = (l){x)v{x, x) + / (f)'{s)v^{s, x) ds ■ ei. 



Clearly, u G C^'"([— a, 00) x Q) and div m = if div u = 0. Further: 

\u\x,x)\ < (IA2I + lAaDll^ziioo + |Ai| ■ |</)(x)m^(0,£) + (1 - (j){x))\\v^\\L^\ 

< (IA2I + \X3\)\\u^\\l^ + |Ai|(|Ai| + IA2I + IAsDIImIl-, 

x)| < (n|A2| + ^^IAsDIIm^IIloo for i = 2,3. 

Since Ai 1, |A2|,A3 0, n|A2| ^ 1 and n^As — > as n — *• 00, the estimate (12.51) 
holds if only n is sufficiently large, which ends the proof. 

We remark that the norm of the operator E blows up when e — >■ 0. Indeed, one 
cannot have e = in (iii) and keep the norm of E bounded. ■ 

The following elementary fact will be often used: 
Lemma 2.3. For any u E L'^ H C°'°(D) there holds: lima;^±oo \\u{x, ■)\\L°°{n) = 0. 

3. The bound on ||m||loo(£)) 

In this section, given c eM., t E [0, 1], a divergence-free vector field v G C^'"{D) and 
a boundedly supported T G C^'°'{D), we consider the following problem: 

— cUx — uAu + rdv ■ Vm + Vp = rTp in D 
(3.1) div M = in D 

u = on dD and lim ||'u(x, ■)||ci(n) = 0. 
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Theorem 3.1. There exists the unique u G fl C^'"(D) solving ^3. 1\) with some 
p e Hf^^{D). Moreover: 

(i) when = then u satisfies: \\u\\l-^[d) < VT||2,2(i^^), 

(ii) when d = 1 then we have: 



\u\\l^(d) < ^^^= yp\Cpw + (^J \p - (0,a;)p^ j \\v\\]^^^^^^\\\/T\\l2(^d)+C\\\/T\\l2 



Above, Cp and Cpw denote, respectively, the Poincare and the Poincare-Wirtinger 
constants offl, while the constant C is independent ofc,T,a,T or v. 

Proof. 1. The hound by VT. Define the quantity: 

1/2 

|p-(0,x)p 



and consider the following vector field with boundedly supported gradient: 

x)=p-ei [ i f{s, ■)ds + p- (0, x) i f{x, ■)• 
JO Jn Jn 

By an easy calculation we see that: 

^\ h f(T .\ 1 n- n. m ^\ X 

dx 



T{x, x)p — Vg(x, x) = ^T(x, ^) ~ J' T{x, p — p- (0, x) -/ -^T{x, ■)ei 



and therefore: 

(3.2) \\fp- VqWmD) < {\p\Cpw + L) II Vf |U2p). 

Recall that the Poincare-Wirtinger constant Cpw on VL is the inverse of the first nonzero 
eigenvalue of the related Neumann problem. By a mollification argument, we may also 
assume that q G C^iD) and that is still satisfied. 

2. Existence of a weak solution. Following the Galerkin method, define: 

V = clmiD) {u e C^{D,R^), div u = O} . 

Clearly, is a Hilbert space with the scalar product {u,w)v = J^^u : Vw. The 

1/2 y 

norms ||u||y := {u,u)y and ||u||//i(i:)) are equivalent in V, in virtue of the Poincare 
inequality in f2, which yields: ||m||l2(d) < Cp||^||y- 

Since is a subspace of H^{D), it is also separable and hence it admits a Hilbert 
(orthonormal) basis {ipn}'^=i G C^(-D) fl V. For each n, let Vn = span {-^i . . . ipn} and 
let Pn : Vn — > Vn be given by: 

Pn{u) = ^ iJi)v-c J u^^Ji + j {v ■ Vu)ipi - T J {T p - Vg)V^i| %IJi. 



8 



MARTA LEWICKA AND PIOTR B. MUCHA 



The operator P„ is continuous and it satisfies: 

{Pn{u),u)v = v\\u\\y — c I UxU + rd i (v-Vu)u — T i {Tp — Vq)u 

Jd Jd Jd 

2Cp 

> z/||M||y - Cp\\Tp- Vq\\L2(D)\\u\\v > when ||M||y = ——\\Tp- Vg||L2(^), 

where we used 2 Jj^u^u = Jjj{\u\^)x = and the nulhty of the trihnear term. 

By Lemma 2.1.4 in p6], there exists G Ki, bounded in V by the above quantity 
and solving: Pniun) = 0. Since V is reflexive, it follows that {«„} converges weakly 
(up to a subsequence) to some u E V such that: 



(3.3) Vw E V u Vu : Vw — c UxW + rd {v ■ \/u)w — t Tpw = 0. 

Jd Jd Jd Jd 

This identity follows first for w = ipn, and then by the density of the linear combinations 
of {V'n} in V- Taking w = u and using (13. 2p we obtain: 

(3.4) \\u\\v < ^ i\p\Cpw + L) ||VT|U2(^). 



3. Regularity. Recall that by de Rham's theorem (see, for example. Proposition 1.1.1 
in [IE]), a necessary and sufficient condition for a distribution field v eT)' that v = Vp 
for some p E V is that {v, w)v = for all w E V. Hence, (13. 3 p implies the first equality 
in (13. ip in the weak sense. By the standard regularity theory [U [16] and in view of 
(13. 4p we may deduce now that the same equality holds in the classical sense and that 
u E H%D), Vp E H\D) (since VT E L'^{D)). Therefore u E C^'^^iD) (for a < 1/4) 
and the boundary conditions in (13. ip follow, together with the asymptotic conditions 
as |a;| — ^ oo, in view of Lemma [2.3[ Next, recalling that T E C^''^{D), the potential 
theory [8] employed to the localized problem and the classical Schauder estimates give 
that u E C^^'^Id). 

4- The bound on cUx- Since Vg has a bounded support and is C^, thus V(p — g) E 
H^{D). Consequently: 

/ UxV{p-q)= I {uV{p-q))x- / div (m(p - g)^) = 0, 
Jd Jd Jd 

because in view of n G H^{D) and V(p — g) G H^{D) one has: 

lim I / |mV(p - g)|(x, ■) + / \u^{p - q)x\{xr)] =0. 
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Integrating the first equality in (13. ip against cu^ on D we obtain: 

11^1*3; 1 1x2(1)) = 1^0 Vm : \/Ux — / cUxirTp — rVg) + rd cUx{v ■ Vu) 
Jd Jd Jd 

< \\cUx\\l^d) (\\Tp- Vg||L2(£,) + d\\v ■ Vu\\L2(^D)j , 
where we have once more used that u G H^{D). Therefore, by (13. 2p : 
(3.5) ||cM:c||l2(£,) < {\p\Cpw + L) \\Vf\\L2(^D) + d\\v ■ Vu\\l2(^d)- 

5. The hound and uniqueness for d = 0. It is now easy to conclude the proof, when 
d = 0. Using the standard elliptic estimates for the Stokes system (12. ip . following from 
the theory in [T], and the Sobolev interpolation inequality, it follows that: 



(D), 



in virtue of (13. 2p . (13. 4p and (13.50 . The constant C is uniform and depends only on 
the geometry of fi, and the constants \p\ and v. Uniqueness of u also follows from the 
above bound. 

6. The case of d = 1. Denote: g = cu^ — rv ■ Vu + {rTp — rVg). By (13.20 and (13.50 
we obtain that: 

\\g\\L^{D) < 2||^||l-'(d)||Vm||l2(£,) + 2 {\p\Cpw + L) ||VT||i2p). 
Therefore, Theorem 12.11 implies: 
2 



\u\\l^(^D) <^=\\v\\^11(^o)\\^u\\l2(^d) 

2 



+ ^ {\p\Cpw + Lf^ II VTlli^/' II VMllii' + Cell Vm|U2p) 



(_D)II ^ "'\\L'2{D) 



which by (13. 4p establishes the result. 



4. The uniform bounds and existence of traveling waves 

In this section we prove the uniform bounds on solutions to system (11.30 . and then 
establish existence of a traveling wave in (11.10 by a Leray-Schauder degree argument. 

Given c G M, r G [0, 1], a divergence-free vector field v G C^'"(i?a) and Z G C^'"{Ra) 
consider first the reaction-advection-diffusion problem: 

- cTx -AT + TvVT = Tf{Z) in 

(4.1) T{-a,x) = 1, T{a,x) = for xefl 

dT ^ 

— -^■(x, x) = for X G [—a, a] and x G dU, 



IV 



(v) 



10 MARTA LEWICKA AND PIOTR B. MUCHA 

together with the following normalization condition, whose eventual role is to single 
out a correct approximation of the traveling wave in T, in the moving frame which 
chooses to have /(T(x, ■)) = for x > 0: 

(4.2) max{T(x,x); x G [0,a], x G 11} = Q^. 

We now recall the bounds on solutions to the above problems, proved in [3] and used 
in [5l[T0]. The right hand sides of (iii), (iv) and (v) follow by re-examining the proofs. 

Theorem 4.1. Let T = Z e C^'''{Ra) satisfy gJ^ and (g^j. Then one has: 

(i) T(x, x) G [0, 1] for all {x, x) G Ra, 

(ii) T{x, x) < 6q for all x > 0,x E Q, 

(iii) |c| < ||^^||L-(i?„) + 2||/'||^{^([o^i]), 

||VT||i2(^^) < 1^1 Q||c - v^\\loo(^r^) + , 
^ fiT)<\n\ (^4||c-t;i|Uoo(«„) + i^. 

Given T G C^'"(i?a) satisfying boundary conditions as in (14. ip . we will consider its 
boundedly supported C^'"{D) extension: 

{T{x, x) for X G [—a, a] 

(l){\x\ — a) ■ (2 — T(—2a — x, x)) for x < —a 
—(f){\x\ — a) ■ T{2a — x, x) for x > a. 

Here (f) G C~(R, [0, 1]) satisfies 0(x) = 1 for x < 1/3, 0(x) = for x > 2/3 and 
<4. 

Also, for a divergence-free v G C^'°'{Ra), let f G C^'°(D) be its divergence-free 
extension, given in Theorem 12. 2[ 

Theorem 4.2. LetceR, t e [0, 1], T = Z G Ci'°(i?a), v G Ci'°(i?a) and u G Ci'"(D) 
satisfy ^3. (JjJ^j H-^ j '"'^^^ ^ ^^'^^^ ^ defined as above. Moreover, let u\r^ = v and 
assume that either d = 0, or d = 1 and holds. Then, for large a: 

|c| + \\T\\ch<^(R,) + ||M||c2..(ij„) + ||VT||l2(^,^) + \\u\\h3(r,) + / /(T) < C, 

jRa 

where C is a numeric constant independent on a, t and the estimated quantities. 
Proof. By Theorem 14.11 (iii) and (iv) we obtain: 

IIVTII^^^^^^ < 7|l]|(||n||ioo(B) + ||/'||i{^([o^i])) + 
On the other hand, the boundary conditions for T imply that, for large a: 
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Consequently, for = the uniform bound on ||m||loo(d) follows by Theorem 13.11 (i). 

When d = 1 then Theorem 13.11 (ii) yields the same bound, under condition (11.41) . 
Indeed, let e > be such that the quantity in the left hand side of (11. 4p is strictly 
smaller than 1 / y/1 + e. Then: 

hh^iD) < VTT7^^ (^\p\Cpw + (^j- \p- (o,x)|2^ j X 

(4.4) X {Vu\Q\'^'{\\u\C^^^ + \\f'\Cao,w) + C) \\u\\'/l 



L°°{D) 



< 9||m|U-(d) + + l)(||/'||/oo([o,i]) + 1), 

1 1/2 

Il-{[o,i]) 



for some q G (0,1). Hence ||m||j;^oo(£)) < C(||/'||^{^(ro id + l) and so, recalling (13.41) and 



Theorem 13. H we have: 
(4.5) |c| + Ml^^d) + Ml.^D) + / f{T) < C(||/'||i/J([o,i]) + 1). 

J Ra 

In (14. 4p and (14.51) the constant C is independent of a, r, the nonlinearity / and the 
estimated quantities. 

The uniform bounds on ||u||/^2(£)), |c|, || VT||i2(/j^) and f(T) follow by Theorem 
14. H (13.41) and (12.41) . Now, the standard local elliptic estimates for the Stokes system 
(12. ip (see yj, also compare flO\) imply that: 



\u 



\m{D) < C{\\Vg\\L2(D) + \\u\\m{D)) 



Taking g = cux — rv ■ Vn + rTp, we obtain the uniform bound on ||m||/^3(£,). The 
bounds on ||m||c2,q(£)) and IjT'llci.aiiJa) follow by Holder's estimates for system (13.11) and 
flirjl P E]. The proof is done. 

We remark that our result does not imply smallness of C in the uniform bound. In 
particular, C depends on the constant C^^n from Theorem [TT] , which can be arbitrarily 
large. ■ 



We finally have: 

Proof of Theorem 11.11 For every sufficiently large a > 0, consider an operator: 

Ka-.Rx Ci'"(i?J X C]'"(i?J X [0, 1] ^ M X C'^'iRa) X C]'"(i?J, 

where C^°'{Ra) stands for the Banach space of the divergence- free, C^'" regular vector 
fields on the compact domain Ra. Define: 

Ka{c,Z,v,T) := (c - 6^0 + max{T(x,x); xG[0,a], x eVt},T,u\R^), 
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where T is the solution to (14.11) . and u solves (13. ip are known also with T and v defined 
as in (14.31) and Theorem 12.21 

The operator Ka is continuous, compact fl] and all its fixed points (c, T, v) such that 
Ka{c, T, V, t) = (c, T, v) for some r G [0, 1] are uniformly bounded, in view of Theorem 
14.21 We may now employ the Leray-Schauder degree theory, as in [31 [5l [10], to obtain 
the existence of a fixed point of Ka{-, ■, ■, 1), since the degree of the map Ka{-, ■, ■, 0) is 
nonzero. This fixed point (c",T",f") again satisfies the bounds in Theorem 14.21 

By a bootstrap argument we moreover obtain the uniform bound on ||T°||c2,a(i{^_^). 
One may thus choose a sequence a„ — oo such that c„ := c""" converges to some 
c e M, and T„ := T"", f„ := v"" converge in Cfj°(D) to some T,u E C^'^'iD). Further, 
u G H'^{D) n C^'"(D) and hence the first convergence in (11.51) follows. Since VT G 
n C'''"(D), we obtain the other convergence in view of Lemma [2.31 

The posit ivity of the propagation speed c and the existence of the right and left 
limits of T, together with the statement in (ii) follow exactly as in [10]. ■ 



The following lemma improves on the result in |10j, where a sufficient condition for 
the left limit 6'_ of T to be 1 required a cubic bound: f(T) < k[(T — Oq)^]^. 

Lemma 5.1. In the setting of Theorem \l.l[ if moreover the nonlinearity satisfies: 



(5.1) /(T) < k[{T - 9o)+]' and k (l + ||/|Uoc([o,i]) + ||/'|U-([o,i])) < Cn, 



then 6_ = 1. Here > is a constant, depending only on u, p and fi. 

Proof. 1. In the course of the proof, C will denote any positive constant depending 
only on p and Vl. Integrating the temperature equation in (11.31) against T and A on 
D yields, respectively: 



where in (15.31) we used that f^T^AT = 0. The interpolation. Holder and Sobolev 
inequalities imply that: 



5. A SUFFICIENT CONDITION FOR 6 



1 



(5.2) 




(5.3) 



\\AT\\l2^d) < ||/(T)|U2(z,) + ||n-VT|U.(^), 




Now, taking e above sufficiently small and introducing (15.31) and (15. 4p into: 

||VT||^,i(^) < C {\\VT\\l2^d) + ||Ar|U2(^)) 
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we obtain, in view of (15. 2p : 

\\^T\\HiiD) < C {\\VT\\mD) + \\f{T)\\L-i^D) + lkl|^i(D)l|Vr|U2(^) 

<c{l + \\f\\'il^^,^,^) + M?miD)) [IJ(T)) ' ■ 
By (14.51) and convergences established in the proof of Theorem II. ![ this imphes: 

(5.5) II vr||^.(,,) < c (i + \\f\\];i + lift/') /(T)) . 

2. Now, for every a; G M denote M(x) = max^jg^ ^'(x, x), m{x) = min^gQ T(x, x) 
and notice that m{x) is non increasing. This can be proved for each T„ on using 
the maximum principle. Passing with n to oo, one obtains the same result in the limit. 

We now argue by contradiction. If 6^ < Oq then m{x) < 9q for every x G M and: 

/ [(T - 9a)+f < \n\ / |M(x) - m(x)|' dx 

J D J -oo 

/+00 /• 
||T(x, T(x, ■)\\L^(^n) dx < C\\VT\\l,^oy 

■oo Jo. 

Together with (15.51) and the assumption in (15.11) the above yields: 

/ [{T - e,)^f <c{i + ii/iUoo + ii/'iu.) k ![{T- e,)+]\ 

J D J D 

which by the assumption on k implies that both sides above must be zero. Conse- 
quently, be /(T) = and one can deduce (see [31 [10]) that T = as well, contradicting 
the results of Theorem 1 1.1[ ■ 

The condition (15.11) seems to be artificial and we believe that it can be further relaxed 
or even omitted altogether, for the wave (T, u) obtained in the limiting procedure of 
Theorem 11.11 
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